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Abstract. We present an analytical and numerical study of the two-dimensional capillary-driven thin film 
equation. In particular, we focus on the asymptotics of its solutions. Linearising the equation enables us to derive 
the associated Green's function and therefore obtain a complete set of solutions. Moreover, we show that the 
solution for any summable initial profile uniformly converges in time towards a universal self-similar attractor that 
is precisely the Green's function multiplied by the initial algebraic volume. Finally, a numerical study enables us 
to conjecture the extension of our results to the nonlinear equation. 
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1. Introduction. Dimensional analysis is well understood through the Vaschy-Buckingham 
LT theorem |13j . Historically, this framework has already led to remarkably important results such 
as the expression of the hydrodynamical drag force on a sphere by Reynolds |41j , the theory of 
turbulence by Kolmogorov [25]. I26|. and the estimation of the nuclear explosion power by Taylor 
[501 [FT] , Moreover, dimensional analysis is directly connected to the fundamental concepts of 
scaling and self-similarity, that appear in numerous situations such as fractals and diffusion |20| . 

The powerful theory of intermediate asymptotics developped in particular by Barenblatt goes 
one step further in understanding the deeper meaning of self-similarity [4j. In nonlinear problems, 
one may wonder what is the interest of finding exact particular solutions as there is no principle 
of superposition. Nonetheless, in certain cases, the self-similar solutions obtained for idealised 
problems, or idealised initial conditions, represent the intermediate asymptotic regimes of the 
solutions of more general non-idealised problems. Following Zeldovich in the foreword of jl], one 
might even say that intermediate asymptotics is the key that somehow replaces the superposition 
principle in nonlinear physics. In other words, by paying the price of a loss of information at 
short times, one obtains a certain generality at intermediate times. Even so, such an asymptotic 
behaviour must be proved for any given initial condition, which often turns out to be a difficult 
task. 

Furthermore, in the afterword of his book Barenblatt states [4]: "However, there exist many 
problems of recognised importance where this technique has not yet been fully explained, but for 
which results of substantial value can be expected from its application. " One of these open problems 
is precisely the capillary -driven thin film equation of interest [51 1171 [55] : 

d T H + d x (H 3 d x H) = . (1.1) 

It governs the capillary evolution of the profile H(X, T) of the free surface of a thin viscous liquid 
film: as soon as the curvature is nonconstant, this profile is unstable as the Laplace pressure drives 
a flow that is mediated by viscosity. Despite many efforts, this equation remains only partially 
solved. Nevertheless, in the past few years, several analytical |11[ 157] HI] an d numerical studies 
[8j 04] have been performed in order to gain insights into this mathematical problem. 



It should be stressed that Eq. (1.1) is of tremendous importance in a variety of scientific 
fields such as polymer physics, physiology, biophysics, micro-electronics, surface chemistry, ther- 
modynamics and hydrodynamics, since thin films are involved in modern mechanical and optical 
engineering processes, through lubrication, paints and coating. Gaining a complete understand- 
ing of these systems is a key step towards the development of molecular electronics, biomimctics, 
superadhesion and self-cleaning surfaces. Moreover, this equation may be crucial for understand- 
ing the nanorheology of ultra-thin polymer films, for which enhanced mobility effects have been 
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predicted [12], before being related to entanglements networks [48], and observed in various ex- 
perimental configurations such as: confinement (31 EHl US HZ] > glassy state [TU], and dewetting 
onto slippery substrates |51|3B]. It may also govern the surface instabilities and pattern formation 
[2, 15, 33 , through the film preparation by spin- coating [6 39, 40, 49 . From all these examples, we 
understand the necessity of exploring further the solutions of the capillary-driven thin film equa- 
tion. In particular, as for the diffusion equation [20J, it would be interesting to characterise the 
convergence of the solutions to some asymptotic self-similar regimes [13], as well as the existence 
of possible exotic self-similarities [46] , 

The present article is divided into three parts. In the first one, we recall the main ingredients 
of the physical model that describes two-dimensional capillary-driven flows. In the second one, 
we linearise the governing equation and we derive the Green's function that enables to calculate 
the general solution of the linear problem for any summablc initial condition. In particular, we 
study the self-similar asymptotics of this solution. In the third part, we extend these ideas to the 
nonlinear equation through numerical solutions. 



2. Physical model. In this part, we describe the model and derive the capillary-driven 
thin film equation, within the lubrication approximation, for two-dimensional viscous flows. The 
evolution in time t of a thin viscous film described by a profile of height z = h(x,t) can be 
understood from the Laplace surpressure p(x, t), which arises due to curvature at the free interface. 
Considering small curvature gradients, one can write: 



p(x,t) ~ -~/d*h , 



(2.1) 



where 7 is the liquid-air surface tension. In the lubrication approximation, the Stokes equation 
along the x horizontal direction is given by |30| : 



d x p = Tj d 2 v , 



(2.2) 



where v(x,z,t) is the horizontal velocity and rj is the dynamical viscosity. We assume no slip at 
the solid-liquid interface and no stress at the liquid-air interface, so that: 



v\ z=0 = 
d z v\ z=h = 



(2.3a) 
(2.3b) 



The pressure being independent of the vertical coordinate z, Eq. (2.2 1 together with Eq. (2.3 1 lead 
to a Poiseuille flow of the form: 



v(x, z,t) = -* 1 (z 2 - 2hz) d x p 
Zrj 



Conservation of volume can be expressed as: 

d t h = -d x 



dz v . 



(2.4) 



(2.5) 



Combining Eqs. (2.1 1, (2.4 1 and (2.5 1 leads to the two-dimensional capillary-driven thin film 
equation: 



d t h 



2_ 

'■iq 



d x (h 3 d*ti) = 



Finally, Eq. (2.6 1 can be nondimensionalised by letting: 

h = H h 
x = X ho 

"irjho 



t = T 
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(2.6) 



(2.7a) 
(2.7b) 

(2.7c) 
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where ho is a typical length scale of the problem. This leads to the dimensionless equation 
introduced above in Eq. ( |1.1[ ). This equation can be linearised by letting: 

H(X,T) = 1 + A(X,T) , (2.8) 

and by assuming ACl: 

d T A + d^A = 0. (2.9) 
The linear case corresponds to a situation in which the surface of a flat film is only slightly per- 



turbed. Note that Eq. (2.9 1 also describes surface diffusion phenomena leading to flattening of 
solid surfaces |35l 154]. kinetic growth [28, 25] [52], or grooving |34l |4"3"]. Therefore, the following 
results apply to a broader class of physical situations. 



3. Solving the linear equation. In this part, we solve Eq. (2.9) and we explore its solu- 



tions. The linear study is divided into four paragraphs. In the first one, we derive the Green's 
function and show that it is self-similar at all times. In the second one, we give the general formal 
solution and exploit it through two particular canonical examples. In the third one, we study the 
uniform convergence in time of the general solution towards the self-similar Green's function. In 
the fourth one, we conclude the discussion on the linear case with some general remarks. 



3.1. Green's function and self-similarity. Since Eq. ( |2.9| is a linear partial differential 
equation, it can be solved by calculating the Green's function Q(X,T). This object is defined by: 

[dr + d*] G(X,T) = S(X,T) , 

where 5 denotes the Dirac distribution. Fourier transforms are defined as follows: 

1 

W? 
1 

dXdT Q(X, T) e - l(nT+KX) 



G(X,T) = 
G{K,T) = 

G{K,Sl) = 



dKdSig(K,Sl)e^ nT+KX ^ 
dSl G(K, SI) e inT 



(3.1) 

(3.2a) 
(3.2b) 
(3.2c) 



Writing the Fourier transform, as defined in Eq. ( |3.2| ;), of Eq. (3.1 1 leads to: 

1 



G{K, SI) 



iSl + K 4 



Using Eq. (3.2 :>) and Eq. (3.3 I, one obtains 

'0 
Res 



G{K,T) 



HIT 



SI - iK 4 



iK l 



if T < 
if T > 



(3.3) 



(3.4) 



where Res(/;z*) denotes the complex residue of the function / at z — z* . Note that Eq. (3.4 1 



is consistent with causality. Henceforth, we shall restrict to T > in the following. Finally, 



expressing the residue and performing the inverse Fourier transform, defined in Eq. (3.2i), one 
obtains the Green's function: 



G(X,T) = lJdK 



e~ KT e iKX 



(3.5) 



which is consistent with previous studies [111 I34j. Then, at finite time, let us change variables 
through: 



X = UT 1 ^ 
K = QT- 1/4 



(3.6a) 
(3.6b) 



M. BENZAQUEN, T. SALEZ AND E. RAPHAEL 




Fig. 3.1. Auxilliary function <f>(U) = T 1 / 4 Q(U,T) given in Eq. | |3.9| , where Q(U,T) is the Green's function 
of the dimensionless linear two-dimensional capillary- driven thin film equation given in Eq. (12. 9| . 



This gives: 



T 

where we introduced the auxilliary function: 



G{U,T) = 1 ± Ii <j>{U), (3.7) 



<t>(U) = ~ / dQe- Q4 e lQU . (3.: 

2-7T J 



The Green's function is thus self-similar of the first kind [3], at all times. Furthermore, the 
function <fi is given for all U G K by: 

^>=^(!) »-({i!}.(?) 4 )-^-^(!) -<{!.|}.(?) 4 ) 

where the (0, 2)-hypergeometric function is defined as [1, 21 : 

o^({M},-) = £^^, (3.10) 



with the Pochhammcr notation (.)t for the rising factorial. The function <f> is plotted in Fig. 3.1 



Apart from the oscillatory behaviour which is directly related to the fourth spatial derivative of 



Eq. (2.9 1, this solution is qualitatively close to the point-source solution of the heat equation [20 , 
for which the same analytical treatment would lead to the well-known Green's function and to its 
specific self-similar variable XT~ X I 2 . 



M. BENZAQUEN, T. SALEZ AND E. RAPHAEL 



5 




Fig. 3.2. Normalised analytical solution (see Eq. | |3.12| ) of the diraensionless linear two-dimensional 
capillary- driven thin film equation (see Eq. ( |2,9| l ) as a function of the self-similar variable U introduced in Eq. ( |3.6| . 
The initial profile is defined as a gate function of volume unity, as shown in the inset. The solution is plotted for 
different dimensionless times and compared to its normalised asymptotic attractor given in Eq. (|3.9[l . 



3.2. General solution. For a given summabl^j] initial condition, A(X,0) = Ao(X), the 
solution A(X, T) of Eq. (2.9 1 is given by the spatial convolution of Aq(X) and the Green's function: 



A(X,T) = J dYG(X -Y,T) AoOO . 



(3.11) 



Interestingly, Eq. (3.11) implies that the Green's function of the problem is exactly the point- 
source solution obtained from an initial Dirac spatial distribution: Aq(Y) = S(Y). Let us once 



again change variables through Eq. (3.6l. Then, Eq. (3.5) and Eq. ( 3.11[ ) lead to 



A(U,T) = 



1 



2tt TV4 



dQe 



-Q 4 e *QU 



(3.12) 



The particular case of a stepped initial conditio?^ was studied in detail in a previous communi- 
cation |45j . In the present study, we have calculated the solutions for various summable initial 
conditions. We present two of them corresponding to canonical illustrations: a gate function (see 
Fig. 3.2 1 and a gaussian function (see Fig. 3.3 1. In each case, we plot the normalised analytical 
solution of Eq. (2.9 I obtained from Eq. (3.121, as a function of the self-similar variable U intro- 



duced in Eq. (3.6 1, the initial profile at T — being shown in the inset. The solution is plotted for 



different dimensionless times. For comparison, we also plot the normalised function tfi(U)/ (f)(0), as 
given in Eq. (3.9 1. As one can see, both profiles seem to converge in time towards the later function 



1 Along the present article, summable means Lebesgue integrable, which implies in particular that 
m Ao(X) = 0. 

ioo 

2 For comparison, in such a case one gets: 



A(C/,T) = !(l + / d QJ-e-«\ 
2 v J l7T Q 



iQU 



where the dashed integral represents Cauchy's principal value. This solution is self-similar for all T > 0. 



G 



M. BENZAQUEN, T. SALEZ AND E. RAPHAEL 




Fig. 3.3. Normalised analytical solution (see Eq. | |3.12| ) of the dimensionless linear two-dimensional 
capillary- driven thin film equation (see Eq. ( |2.9| ) as a function of the self-similar variable U introduced in Eq. | |3.6| . 
The initial profile is defined as a gaussian function of volume unity, as shown in the inset. The solution is plotted 
for different dimensionless times and compared to its normalised asymptotic attractor given in Eq. \3.9\. 



that we shall henceforward call a universal self-similar attractor. This convergence statement will 
be addressed in the next paragraph. 



3.3. Uniform convergence to the self-similar attractor. For clarity, we shall restrict 
here to summable initial profiles that have non-zero algebraic volume, that is: 



Mo = J dXA (X) 



(3.13) 



The extension of the following results to the specific case of zero initial algebraic volume is under- 
stood as well and will be addressed in the next paragraph. For the time being, let us introduce 
the function: 



T i/4 

MT)= — A(U,T) 



For all T > 0, and for all U G M, one has 

1 



\f(U,T) - 4>(U)\ < 



< 



2tt\M \ 
1 

2n\M \ 



dQ 



e -Q 4 e iQU 



dY [e-^ Y l Tl/l - l) A (Y) 



dQe~ 



dY 



-iQY/T 1 / 4 



\MY)\ 



By introducing: 

a(T) 

one gets for all T > 0: 



2tt|X c 



dQe~ 



dY 



|Ao(F)| 



\\MT) -<j>{U) |L„<o(T) . 



(3.14) 

(3.15) 
(3.16) 

(3.17) 
(3.18) 
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where ||...||oo,c/ is the uniform nornxjwith respect to the variable U. In order to conclude that the 
function / is uniformly convergentrfin time towards <j>, it remains to show that: lim a(T) = 0. 

For this purpose, let us consider the auxilliary function h(Y, T) defined by: 

h(Y,T)=(e-^ Tl/4 -l) A (Y), (3.19) 

for all Y £ K, and for all T > 0. This function naturally converges to the zero function when 
T — » oo . In addition, for all T > 0, and for all Y £ K, one has: 

\h(Y,T)\ < h(Y) , (3.20) 

where we introduced the summablc function: 

h(Y) = 2 \A (Y)\ . (3.21) 

Invoking Lebesgue's continuity theorem leads to: 

lim / dYh(Y,T) = . (3.22) 

Similarly, we then consider the second auxilliary function g defined by: 

g(Q,T) = e~ Ql J dYh(Y,T) , (3.23) 

for all Q £ M, and for all T > 0. It converges to the zero function when T — > oo, according to 
Eq. ( |3.22| . In addition, using Eq. fl3.20| and Eq. ( |3.21| , for all T > 0, and for all Q £ K, one has: 



\g(Q,T)\ <g(Q) , (3.24) 
where we introduced another summable function: 

g(Q) =2e-^ J dY \A (Y)\ . (3.25) 

Once again, invoking Lebesgue's continuity theorem leads to: 

lim a(T) = . (3.26) 



In summary, we thus demonstrated that the solution A(U,T) of Eq. (2.9), given in Eq. (3.12) for 
any summable initial profile, always converges uniformly in time towards the universal self-similar 
attractor defined as: 

A 0o (U,T) = ^d>(U)- (3-27) 



We refer one more time to Fig. |3.2| and Fig. |3.3| for illustration of this result with a gate function 
and a gaussian function as initial profiles, respectively. The crucial point in these graphs is that 
the attractor is identical for the two summable initial profiles. 



3 |[/|[cx>, :c =sup{|/(a;,2/)|,a ; eR}. 

4 In the present article, a function of two variables f(x, y) is said to be uniformly convergent with respect to y 
if lim H/lloo,* =0. 

y — >oo 
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3.4. Remarks. In order to conclude the discussion on the linear case, we enumerate five 
general remarks below. 

• The latter calculations may be extended to all equations of the form: 



[di 



9l m l A{X,T) 







(3.28) 



with m e N*. Indeed, apart from the well known heat equation [20], one can obtain the 
Green's function for higher even orders of the spatial derivative and extract analogous 
conclusions. However, the odd orders being free from dissipation are expected to lead to 
fundamentally different mathematical solutions. 

According to our primary interest in thin films, we proved the uniform convergence of 
any summable solution towards the self-similar attractor by explicitly writing the Green's 
function of Eq. ( |2.9[ ). However, this result is more general. In fact, let us consider any 
diffusive-like linear partial differential equation of two variables for which the Green's 
function G(X, T) is self-similar of the form: 



G(X,T) 



1 



X 



(3.29) 



where a and /3 are strictly positive real numbers, and where <f) is a bounded function on 



Then, the uniform convergence is straightforward to demonstrate, as soon as Eq. (3.111 



and Eq. (3.13) remain satisfied, with a summable initial profile. 

Another interesting feature is the following. Let us call A[a ](A, T) the solution of 
Eq. (2.9 1, for an initial profile A (A). Then, the evolution A^^(X, T) of the first deriva- 
tive Aq of the previous initial condition is simply given by SxA[a ] (X, T). In terms of 
the self-similar variable introduced in Eq. (3.6), one gets: 



A[Ai](C/,T) 



T l/4 



d v A [Ao] (U,T) 



(3.30) 



This notably makes the link with our previous work on the particular case of a stepped 
initial condition |45| . Indeed, the Dirac distribution being the first derivative of the 
Heaviside's distribution, the Green's function is thus simply given by the derivative of the 
solution obtained for a stepped initial condition. Generalised to higher order derivatives, 
this relation naturally becomes: 



1 



jn/4 



A[A n] (£/,T) 



(3.31) 



One may as well wonder what happens in the particular case of zero algebraic volume^] 
that is when Mq = 0. The answer is that there is still an attractive self-similar regime 
given by the first non-zero derivative cf>^(U), under the condition that this quantity is 
summable. 

At last, we have seen that for any summable initial condition there is long-term self- 



similarity of the general solution of Eq. (2.9 1. The question arrises to know whether 



there exists some specific summable initial conditions that generate solutions that are 
self-similar at all times. Interestingly, the solution obtained from the Heaviside initial 
condition was shown to be self-similar at all times for other boundary limits [45] , In the 
present case, let us impose the following constraint at all times T > 0: 

A(C/,T) = T a/4 F(U) , (3.32) 



5 One could for instance imagine a dip followed by a bump of identical shape, as a pathological initial profile. 
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with U as defined in Eq. (3.6 I. By changing variables in Eq. (3.121, it is straightforward 



to see that the initial profile must be homogeneous of degree a, meaning that for all real 
numbers k and Y, one has: 



A (kY) = k a A (Y) . 
Finally, using Eq. (3.131 and Eq. ( |3.33[ ), the algebraic volume V satisfies: 



V = J dX A(X,T) 

= T (a+i)/4 J dUF(U) , 



(3.33) 

(3.34) 
(3.35) 



which implies that a = — 1, since V = Mq by volume conservation. Therefore, the only 
initial profile with finite algrebraic volume that exhibits the self-similarity of Eq. (3.321 
at all times is the Dirac distribution. Recalling the remark made after Eq. (3.11), this 



means that the Green's function is the only summable solution that is self-similar at all 



times, with the definition of Eq. (3.321 



4. Extension to the nonlinear equation. In this last part, we extend the previous results 
to the nonlinear case through a numerical scheme. Thus, we consider the nonlinear partial dif- 
ferential equation given in Eq. (1.1 1. This equation has not been solved analytically yet, but we 
recently solved it numerically in various geometries [44] , The numerical procedure we used is a 
finite difference method developed in |8|I53|. It ensures capillary energy and entropy dissipation, as 
required from [7]. In addition to volume conservation, it has been shown that this method ensures 
positivity of the height profile H(X, T) [53]. Using this numerical scheme, we verified the existence 
of an attractive self-similar regime for several arbitrary initial profiles, with compact support as 
required from the algorithm. For instance, we plot in Fig. |4.l|the normalised numerical solution 




Fig. 4.1. Normalised numerical solution of the dimensionless nonlinear two-dimensional thin film equation 
given in Eq. (TlTJ, as a function of the self-similar variable U defined in Eq. ( |3 ■ 6| , for an initial profile given by an 
arbitrary function with compact support, as shown in the inset. The solution is plotted for different dimensionless 
times and compared to the normalised universal asymptotic attractor of the linear case given in Eq. ||3.9| . 
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of Eq. (fTTTJ function of the self-similar variable U denned in Eq. ( |3.6[ ) , the initial profile being 
shown in the inset. The solution is plotted for different dimensionless times. We also plot for 
comparison the normalised function 4>(U)/4>(0), as given in Eq. (3.9 K What is remarkable here 
is that the attractor of the nonlinear case is actually equal to the function <p(U) / <p(0) obtained 
analytically in the linear case. Even though this may seem unexpected, it turns out to be quite 
natural when one realises that as time goes, the initial profile progressively collapses towards the 
flat film equilibrium shape, thus bringing Eq. (1.1) closer to its linearised form of Eq. (2.9). 



5. Conclusion. We reported on the asymptotics of the two-dimensional capillary-driven thin 
film equation for an arbitrary summable initial profile. First, we derived an analytical solution 
of the linearised equation. The solution was obtained by seeking the Green's function, which was 
found to be given by a combination of generalised hypergeometric functions. We then proved that 
any summable initial condition lead to the uniform convergence in time of the solution towards 
a universal self-similar attractor that is proportional to the Green's function of the problem; the 
proportionality factor being equal to the initial algebraic volume. This result appears to be a more 
general result in diffusive-like processes that are characterised by a self-similar Green's function 
with decaying amplitude. At last, we were able to conjecture from compact-support numerical 
results, as well as to justify on a physical basis, the extension of this convergence behaviour to 
the nonlinear equation. The important outcome is that the universal self-similar attractor of the 
nonlinear case appears to be precisely the one of the linear case. 

The recent excellent agreements between thin film theories and experiments with stepped 
polymer films [3"Tl I3"2l FH1 H5] , and polymer droplets on identical films |16l HI] , are very encour- 
aging for the physical relevance of the present analysis. Thus, experimental implementation with 
viscous nanofilms should be performed in near future. This theoretical work may also be extended 
to other thin film equations [H] , such as the one describing flows in Hele-Shaw cells [I] and the one 
governing gravity-driven flows |18| 123] [27] . One could also address coalescence phenomena |22l [42] 
and various diffusive processes, in an identical way. Moreover, it should be feasible to exhibit the 
necessary conditions acting on a given partial differential equation for showing such a uniform 
convergence behaviour. Finally, characterising further the conditions for finite-time convergence 
may be of fundamental importance for the field of intermediate asymptotics, and thus for nonlin- 
ear physics in general. 
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